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Abstract 

We present a novel flow instability that can arise in thin films of cytoskeletal fluids if the friction 
with the substrate on which the film lies is sufficiently strong. We consider a two dimensional, 
membrane-bound fragment containing actin filaments that is perturbed from its initially circular 
state, where actin polymerizes at the edge and flows radially inward while depolymerizing in the 
fragment. Performing a linear stability analysis of the initial state due to perturbations of the 
fragment boundary, we find, in the limit of very large friction, that the perturbed actin velocity 
and pressure fields obey the very same laws governing the viscous fingering instability of an interface 
between immiscible fluids in a Hele-Shaw cell. A feature of this instability that is remarkable in 
the context of cell motility, is that its existence is independent of the strength of the interaction 
between cytoskeletal filaments and myosin motors, and moreover that it is completely driven by 
the free energy of actin polymerization at the fragment edge. 

PACS numbers: 87.17.Jj, 87.17.Rt, 61.30.-v 
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Directed motion and shape change allow cells to respond to their environment and play 
central roles in many biological processes such as embryonic development, wound healing, 
and formation of cancer metastases. Almost universally, crawling of cells on a surface or 
extracellular matrix involves the protrusion of a thin leading edge, the lamellipodium, driven 
by the polymerization of the actin cytoskeleton, the adhesion to the substrate via specific 
proteins and molecular motor-enabled contraction of the cytoskeleton to translocate the 
traUing cell body fl. Remarkably, physical units far simpler than euUaryofc cells can 
display self- sustained motion: the work of Ref. [2j shows that nearly flat cell fragments, 
containing only actin cytoskeleton and myosin II motors enclosed by a plasma membrane, 
are able to perform polymerization/contraction-driven motion. Furthermore the fragments 
can spontaneously switch between motile and non-motile states. The observations of Ref. ^ 
have led us to study theoretically actin-driven motility and shape dynamics in these simpler 
systems with few structural elements and few measurable parameters. 

The actin cytoskeleton is a highly complex medium: it is polar as actin polymerizes at 
its "plus" end, facing the membrane abutting the lamellipodium; it is viscoelastic; and it 
is active and driven out of equilibrium by ATP hydrolysis, needed for continuous polymer- 
ization (treadmilling) and to generate myosin motor-induced stresses. Recently, a generic 
hydrodynamic theory has been developed to describe active, polar media |3|. Using a small 
number of phenomenological parameters, it can account for a number of motility phenomena 
due to the coupling of actin filaments to myosin activity 0, . The interaction of the cell 
with its environment is also very important in describing motility; for example, it has been 
demonstrated that cells crawling on a heterogeneous substrate tend to migrate to regions 
of greater substrate adhesion [Tj] and and greater substrate rigidity [8]. Cytoskeletal actin 
in a cell or cell fragment is able to transmit forces to its substrate through transmembrane 
proter. MU, namely .nte.nn. wMch Mnd reve^My to the .b.trate. In general, .n- 
tegrins cluster to form focal adhesions whose size and mechanical properties are determined 
by chemical and mechanical cues and can regulate the force that they exert. If, however, the 
actin velocity relative to the substrate is small compared to a/r, where a is a molecular size 
and r is the average time during which an integrin remains bound, then the force exerted 
by the moving filaments on the substrate can be expressed as a friction force, proportional 

to the actin velocity fly 

In this Letter, we demonstrate that polymerization and large friction forces are sufficient 
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to destabilize an initially stationary, circular cell fragment. We start by considering a very 
simplified model of actin cytoskeletal flow in a cell fragment, as shown schematically in Fig.[Tl 
The fragment is very thin and there is no flow or any spatial dependence in the 2;-direction; 
also we consider that the thickness of the fragment is constant. We further assume that the 
material in the fragment can be modeled as a single fluid: this implies that we treat only 
the flow of actin, and assume that the cytosol (including free actin monomers) is stationary 
relative to the substrate. Furthermore, we treat the cytoskeleton as an incompressible liquid, 
ignoring the elastic response that occurs on times shorter than the viscoelastic relaxation 
time. 

Actin polymerization is regulated by proteins such as those of the Wiskott-Aldrich syn- 
drome family (WASP), which localize in the cell membrane at the fragment edge jl4]. For 
the purposes of this work, it is sufficient to assume that actin polymerizes only at the frag- 
ment edge and in the direction normal to the boundary. Newly polymerized actin flows away 
from the fragment edge by treadmilling, due the turnover of free actin monomers back to 
the edge for further polymerization that is enabled by actin depolymerization in the bulk. 
For simplicity, we assume that the filament depolymerization is spatially uniform and occurs 
at a rate proportional to the filament density. 

These simplifications imply that the actin fiow in the unperturbed, circular state induced 
by localized polymerization and uniform depolymerization is imposed by the continuity equa- 
tion V • (pvo) = —kdp, where p is the actin filament density and kd is the depolymerization 
rate. The assumption of incompressibility directly leads to the radially-directed treadmilling 
speed 

^0 = -yr. (1) 

Note that in the stationary state continuity requires that ^-^o = Vp, where Rq is the 
unperturbed fragment radius and Vp is the polymerization velocity. 

The cytoskeleton dynamics is described by the hydrodynamic equations for active polar 
gels of Ref . ^ , which themselves are a generalization of the hydrodynamics of liquid crys- 
tals [15I, [1^ modified to account for the coupling between stresses and active motors as well 
as actin polarization and motors. We can, however, proceed by considering that the friction 
with the substrate only couples directly to the actin fiow and not to the polarization. The 
viscous fingering instability that will be seen shortly to be driven by edge polymerization 
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FIG. 1: (a) Schematic drawing of actin cytoskeleton in the unperturbed, radial state. The direction 
of average actin polarization is radial, and is indicated by p = f . (b) A side view of the fragment. 

and bulk depolymerization in the large friction limit can then be most easily illustrated 
by neglecting the dynamics of the polarization field, p, and assuming this quantity is fixed 
along the radial direction so that p = r throughout. 

Ignoring the dynamics of the polarization field, the constitutive laws of Ref. Sj reduce 
to those for an isotropic, viscous fluid of viscosity rj, augmented by an active term in the 
deviatory stress component a^r that reflects the myosin-mediated interaction between actin 
filaments that are nearly aligned; that is, a = 2r]u — CA/irr, where i] is the actin viscosity; 
u is the velocity gradient tensor; A/i is the chemical potential difference between ATP 
and its hydrolysis products; and where for contractile motors the activity coefficient ( is 
negative [17]. The constitutive laws are completed at low Reynolds numbers by the force 
balance V-a = VP+^v, where ^ is the friction coefficient between the cytoskeletal filaments 
and the substrate. Scaling lengths by Rq, times by 1/kd, and stresses and pressures by rjkd 
(keeping the same variable names for the new, dimensionless quantities), it follows that 
viscosity and friction affect the cytoskeletal dynamics through the dimensionless parameter 

= In the limit of very large friction, that is, A ^ 0, the leading term in the force 
balance is simply 
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cell 



V ~ -A^VP, (2) 

ocity satisfies a two dimensional Darcy's law, as it would for the flow in a Hele-Shaw 
|. Based on the available experimental results we find that the quantity A^ is in fact 



quite small. Taking a value of ^ ~ 10^ Pas/ fj? [l9|, 7] ~ 10'' Pa-s [12, l2Q|, and Rq ^ 10 /xm 
we find A^ ^ 10"^. 

We now perturb the edge of the fragment, so that in terms of the polar angle 6 the 
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fragment edge is now at a position R{6,t) = 1 + SR{9,t). For SR{9,t) ^ 1, we perform 
a linear stability analysis by writing 6R{9,t) = '^'^^i SRm(t) cos {m9), and similarly for 
the two components of the perturbed velocity field, 6vr = Ylm=i^'^'r,m{r,t) cos {m9) and 
= J2m=i ^'^9,m{r, t) sin {m9) , and the pressure field 6P = J2m=i ^Pm{r, t) cos {m9). The 
m = mode is excluded from this discussion since it is trivially stable because the quantity 
of actin in the fragment is fixed. In assuming that the depolymerization rate does not 
change as a result of the perturbation and that the filament density remains unchanged, 
it follows that V ■ (5v = and therefore ^f^.m = ^PP^yi^S -^1- © to the 

perturbed quantities 5\ and 5P, we find that V^(5P = and therefore 6Pm{r) ~ r"* and 
6vr^m{r,t) = —dvo^mijit) = ^m(^)^™^^- The Coefficient A^it) can be found be imposing the 
force free condition at the boundary, namely 

dPo 



'^-^ ar 



0, (3) 



r=l 

leading to Am{t) = '^SRmit). The growth rate of the perturbation modes SRmit) is obtained 
by noting that, to linear order in 6R, 

~ 5Vr,ra{Ro) + ^^rn-J^, (4) 

which, using the expression for Amit), gives d6Rm/dt = UmSRm, where the leading order 
growth rate, in units of k^, is 

^„.~^ + 0(A2). (5) 

Note that the mode m = 1, corresponding to an infinitesimal translation of the circular 
fragment, is marginally stable, as required by translational symmetry. 

The linear dispersion relation, Um, is a common feature to a number of Laplacian growth 
problems, for example the viscous fingerin g in stability that occurs at an interface between 



two immiscible liquids in a Hele-Shaw cell 



21| . The physics of the instability is understood 



as follows. The pressure gradient at the edge is dPo^r = l)/dr > 0, and a perturbation with 
6R < 0, for example, requires a perturbed pressure 6P{r = 1) > to keep the boundary 
force-free, to leading order in A^. An excess pressure at the edge relative to the fragment 
center {6P{r = 0) = for m > 0) drives an inward-directed flow, thus amplifying the initial 
negative perturbation. 

Viscosity, surface tension, and motor activity affect the growth rate, Um, at O(A^). In 
short, one obtains from the radial component of the force balance a fourth order ordinary 



differential equation in r for u 
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+ ^ -u 



(2m2 + l) „ (2m2 + l) , 



+ -v! 
r 



m 



-u 



0, 



(6) 



where the primes indicate differentiation with respect to r. The four boundary conditions 
required are that (5v(r = 0) = 0; that the edge of the perturbed fragment is force-free, namely 
^CTnnl'" = 1) = 0, where the subscript n refers to direction normal to the cell fragment; and 
that, neglecting the viscosity of the fragment's surroundings compared with the viscosity 
of the cytoskeleton, the tangential shear satisfies 6atn{r = 1) = 0, where t refers to the 
direction tangent to the perturbed fragment. 

In the limit A — 0, Eq. ([6]) together with the four boundary is a singular perturbation 



22| . a uniformly convergent 



52(t) + C2(t)e^ 



(7) 



problem. Following the boundary layer techniques of Ref. 
approximation to 6vr,m. on the interval < r < 1, valid to O(A^), is 

where -Bo(''^) = ^^Rmif), is the unperturbed growth rate Amit) calculated above which is 
obtained from the boundary condition 5cr„„(r = 1) = at leading order (0(A^^)); The two 
other constants B2(t) and C2(t) are found, respectively, by solving 5cr„„(r = 1) = and 
^citnir = 1) = at next-to-leading order (0(A°)), giving, by way of Eq. (jl]), a growth rate 
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2m 



+ 0(A3). 



(8) 



Equation ([8]) shows that the stabilizing effect of viscosity is proportional to for large m. 
It can be easily seen that the stabilizing effect of the plasma membrane tension also scales 
as (as does the effect of interface tension in viscous fingering instability in the Hele-Shaw 



cell 



2l[|): including membrane tension, the normal stress at the boundary satisfies a two 



dimensional Laplace law 6a„ 
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6H, where 7 is the membrane tension and 



where the mth mode perturbation in the membrane curvature in the (r, 6) plane (ignoring 
changes in curvature in the 2;-direction) is 6H = [th? — l)5RmCO^{m9). The contribution 
of the membrane tension to A2{t) is — A^;^^^m(m^^l). Taking 7 ~ 10~^ N/m [23| as an 



estimate for the membrane tension and A;^ ^ 1 s ^ 2J] it is clear that the stabilizing effect 
of membrane tension is negligible compared to that of actin viscosity. 



Equation ([8]) further shows that the contractile effect of the motors is to stabihze the 
growth of perturbations, proportional to m^. This result depends strongly on the assump- 
tion, made for simplicity, that the filament polarization in the perturbed fragment remains 
everywhere radial and is not a dynamical quantity in the problem. This assumption is 
questionable, yet it may be valid for small m, where changes in membrane curvature are 
small and hence any membrane- act in filament coupling is unlikely to force a significant re- 
orientation of filaments at the leading edge. In any case, Eq. ([8]) shows that the relative 
contribution of myosins to uJm is proportional to Taking C^l^y. ~ —10'^ Pa [sl, this is of 
order ~ 0.1, and therefore small compared to the viscous contribution at O(A^). 

Diffusion of free actin monomers also limits the perturbation growth; however, we may 
consider for now that the diffusion constant, is such that D ^ k^RQ, so that perturbations 
that are area-preserving at first order in 6R, that is, for m > 1, do not affect the essentially 
spatially uniform monomer density, and hence the polymerization rate, Vp. A more careful 
accounting of the effect of diffusion will be considered in a future publication. 

Finally, it might be experimentally useful to have an estimate of the critical value of 
friction, ^c, for which shape perturbations of a cell fragment become unstable. This critical 
value is defined such that for < "^m < and for > '^m > 0. It is conceivable 
that one could observe the onset of growing shape perturbations by plating cell fragments 
on surfaces of varying degrees of adhesiveness or by culturing fragments from cells that 
have been mutated to weaken or strengthen the binding of integrins to the surface 25|, l26 |. 



Equation ([6]) can be solved numerically for different mode numbers m to find the critical 
value ^ where the growth rate becomes positive as a function of motor strength, \(\Afi/rikd] 
see Fig. [2l 

The numerical estimates of C,c given in Fig. [2] are qualitatively consistent with the value 
obtained by setting Um = in the asymptotic growth rate, Eq. ([8]): lower modes are less 
stable as a function of friction and motor activity has a weak effect on the growth of shape 
perturbations. 

In summary, we have found that large substrate friction and the pressure field created by 
treadmilling in an initially circular cell fragment render it linearly unstable. This instability 
has been analyzed here in the limit 7^ 0, where it shows a close correspondence to 
the classic viscous fingering instability in Hele-Shaw cells. We have also shown by direct 
calculation that the effects of membrane tension, actin viscosity, and contractile motors only 



7 



11 

10 
9 

8 
7 

0.2 0.4 0.6 0.8 1.0 r]kd 

FIG. 2: Critical value of friction, ^c, in units of tj/Rq, versus \C\Afi/r]kd for m = 2 (solid line), 
m = 3 (long dashed line), and m = 4 (short dashed line). Surface tension, a/rik^RQ, is taken to be 
zero. 



affect the flow instability at next-to-leading order in the friction. This instability has the 
potential to be highly relevant to the related biophysical problems of cell shape change and 
cell motility, given that it presents a fundamentally hydrodynamic means for cell dynamics, 
independent of complex biochemical signaling and, significantly, of the presence or absence 



ike to study how the instability presented 
^, in which a circular, stationary fragment 



of molecular motors. In future work we would 
here relates to the fragment experiments of Ref. 
could become anisotropic and motile either spontaneously or due to an external mechanical 
force; the reverse transition was observed in these experiments as well. In the context of 
the work presented here, the metastability of the fragments of Ref. |2|] might be explained 
by their friction with the substrate being just below for the linear instability of the mode 
m = 2. We are at present considering this possibility, by studying a possible finite amplitude 
instability that would couple the motile but shape-preserving mode, m = 1, with the shape 
changing modes, m > 1. 
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